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PREFACE 

This  effort  steamed  from  the  Air  Force  Weapons  Laboratory's  (AFWL)  concern 
to  conduct  high-explosive  testing  under  favorable  atmospheric  conditions* 
thus  insuring  a  minimal  impact  to  off-site  coununlties.  Many  people  partici¬ 
pated  In  this  effort. 

Atmospheric  measurements  were  needed  for  airblast  predictions.  Airblaat 
measurements  were  made  to  protect  the  Government  from  invalid  damage  claims 
and  to  add  to  the  general  airblast  propagation  data  base.  Dr  Robert  E.  Relnke 
(AFWL/NTES)  modified  existing  equipment  to  measure  airblast.  Mr  John  A. 
Leverette  (AFWL/NTES)  performed  the  bulk  of  maintenance,  data  reduction  and 
operation  of  the  equipment  for  field  testing.  For  tests  at  Yuma,  Arizona, 
the  U.S.  Army  Atmospheric  Sciences  Laboratory  at  the  Yuma  Proving  Ground  made 
measurements  of  upper  level  wind  velocity  and  temperature.  For  the  larger 
Yuma  tests,  a  weather  team  from  the  Yuma  Marine  Corps  Air  Station  provided  us 
with  low-level  wind  observations  at  the  test  site.  Finally,  Major  Jon  Kahler's 
(USAF,  Retired,  formerly  of  AFWL/WE)  technical  guidance  throughout  this  project 
was  especially  appreciated. 

I  » 

Metric  units  are  used  throughout  this  report.  Exceptions  are  made  for 
some  meteorological  parameters  which  are  used  in  the  most  commonly  reported ' 
format  as  follows:  Windspeed  in  knots,  atmospheric  pressure  in  millibars, 
and  the  height  of  the  weather  data  in  feet.  Conversion  factors  are  as  follows: 
1  knot  -  0.5148  m/s,  1  millibar  ■  102  pascals,  and  1  kft  *  10^  ft  -  304.8  m. 

All  logarithms  are  to  the  base  ten. 

Any  reference  to  specific  brand  names  is  made  for  identification  only 
and  does  not  imply  endorsement  or  criticism  by  AFWL. 
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INTRODUCTION 

The  Staff  Meteorology  Office  of  the  Air  Force  Ueepone  Laboratory  (AFWL/WE) 
Is  tasked  to  monitor  atmospheric  conditions  and  predict  far- field  airblaat 
propagation  for  AFWL  high-explosive  test  progress.  Under  certain  weather 
conditions,  unacceptable  airblaat  nay  propagate  to  nearby  communities  and 
produce  minor  structural  damage  or  excessive  noise  irritation.  Airblaat 
propagation  predictions*  based  on  the  stats  of  the  local  atmosphere  and  the 
explosive  emounts ,  are  needed  for  the  tests  to  safely  proceed  with  the  assur¬ 
ance  of  no  off-site  damage  and  the  minima  of  coomunity  noise  irritation. 

These  AFWL  high-explosive  tests  involve  conventional  weapons  and  chemi¬ 
cal  datonationo  with  explosive  yields  ranging  from  100  to  over  1*000.000  kg. 
Tests  are  conducted  either  locally  at  Kirtland  Air  Force  Base,  New  Mexico*, 
or  at  a  remote  range  near  Yuma*  Arizona.  The  airblaat  prediction  technique 
for  these  tests  must  be  accurate*  incorporate  changing  weather  conditions, 
and  be  adaptable  for  field  operations  at  remote  sites.  This  report  describes 
the  AFWL/WE  airblaat  prediction  model  which  meets  the  criteria,  and  presents 
the  results  of  both  conventional  weapons  and  chemical  high-explosive  test 
series. 
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BACKGROUND 

For  far-field  alrblast  propagation*  tha  atmosphere  acta  as  a  lana  to 
rai  act  tha  alrblast  vava  althar  upward  or  back  to  tho  ground  surface.  In 
tha  far- field,  a  ahock  wave  from  an  explosion  novas  with  tha  spaed  of  an 
acoustic  vava.  Thus*  tha  dagrao  of  refraction  depends  upon  the  speed  of 
sound  in  layers  of  tha  ataosphera  above  tha  ground  surface.  For  tha  purposes 
of  this  report*  tha  use  of  tha  tarn  "sound  spaed"  will  imply  the  total  sound 
speed  in  a  particular  direction.  At  any  given  level,  the  total  sound  spaed 
along  a  specified  direction  is  approximately  equal  to  the  teaperature  depen¬ 
dent  sound  speed  of  the  air  plus  the  wind  velocity  component  as  given  by  the 
equation: 


VS  -  331  *  [1  +  (T/273)  ] 1/2  -  WS  *  C0S(WD  -  DI)  (-1) 

where 

VS  «  Total  sound  speed  in  the  DI  direction  (m/s) 

Dl  -  Selected  direction  of  interest:  Azimuth 
angle*  clockwise  from  true  north  (360  deg) 
as  viewed  from  the  explosive  source  (deg) 

T  «  Air  temperature  (°C) 

WS  -  Wind  speed  (m/s) 

WD  »  Azimuthal  direction  from  which  wind  is 

blowing;  clockwise  from  true  north  (deg) 


A  sound  speed  vertical  profile  ia  constructed  from  temperature  and  wind 
velocity  measured  at  selected  altitudes.  Tho  magnitude  of  alrblast  propaga¬ 
tion  is  dependent  on  the  shape  of  the  sound  speed  vertical  profile  and  the 
size  of  the  difference  between  the  surface  sound  spaed  and  the  values  aloft. 
If  the  sound  speed  decreases  with  altitude,  the  blast  wave  front  will  be 
refracted  aloft  and  sound  amplitude  will  be  diminished  at  the  surface.  Con¬ 
versely,  when  the  sound  speed  at  soma  height  exceeds  the  surface  value,  a 
portion  of  the  wave  front  will  be  refracted  back  to  tha  surface  with  an 
ensuing  sound  enhancement. 
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Ray  tracing  techniques  ara  frequently  usad  to  compute  atmospheric  eirblest 
refraction.  Several  techniques  are  available  vhich  have  been  used  with  success 
at  other  locations.  Table  1  provides  information  on  a  few  of  these  techniques 
and  models.  These  models  are  unsuitable  for  remote  field  work  as  they  require 
either  extensive  calculations  on  a  mainframe  computer  or  a  highly  subjective 
evaluation  of  the  sound  speed  vertical  profile.  This  report  describes  a  rela¬ 
tively  simple  technique,  the  Blast  Operational  Overpressure  Model  (BOOM) ,  for 
predicting  airblast  propagation  at  remote  sites  without  the  benefit  of  main¬ 
frame  computers.  No  claims  are  made  regarding  the  smelts  of  this  technique 
compared  with  others  in  use.  The  BOOM  was  develop* -1  for  use  in  field  opera¬ 
tions  at  remote  locations  and  it  is  capable  of  producing  cccurate  results 
using  a  limited  memory  portable  computer. 
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TABLE  l.  BLAST  PROPAGATION  PREDICTION 


Model 


Method 


1.  Sound  Intensity 
Prediction  System 
(SIPS) 

(Ref.  I) 


Ray  tracing  to  determine 
focal  points.  Capable  of 
considering  terrain. 
Empirical  sound  intensity 
predictions. 


2.  Inverted  FACT 
Model 
(Ref.  2) 


Deterministic.  Uses 
asymptotic  approximation 
to  prevent  infinite 
intensities  in  the 
regions  of  caustics. 


3.  A  Prediction 
Method  for 
Blast  Focusing 
(Graphical) 
(Refs.  3.4  and 
and  S  describe 
similar  methods) 


Velocity  of  sound  profiles 
are  used  with  graphs  to 
determine  focal  point 
location.  Empirical 
overpressure  equations. 
Subjective  graphical 
analysis. 


4.  Focus:  Comput¬ 
erised  Aid  for 
Making  Sound 
Propagation 
Forecasts 
(Ref.  6) 


Ray  tracing.  Computes 
focus  location  and 
overpressure  amplification 
factor 


TECHNIQUES 

Computational 

requirements 

Mainframe 

computer 


Mainframe 

computer 


Sound  Velocity 
profile  and 
overpressure 
predictions  can 
be  calculated 
with  pocket 
calculator. 

Microcomputer 
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PREDICTION  MODEL 

Richard  Lorens  (R«£.  7)  developed  an  eaplrlcel  overpressure  prediction 
model  for  explosions  of  Mark-82*  500-lb  bombs  st  Blood svorth  Island*  Maryland. 
Overpressure  measurements  were  made  25  km  from  the  detonations  of  air-dropped 
Mark  82  bombs  with  air blast  propagation  over  flat*  marshy  terrain  and  open 
water.  Lorens  Incorporated  the  refractive  effects  of  the  atmosphere  into  a 
single  function*  the  Beta  parameter.  This  parameter  represents  the  atmos¬ 
pheric  condition  which  has  the  greatest  effect  on  airblast  refraction;  i.e.* 
the  maximum  difference  in  the  speed  of  sound  between  the  surface  and  the  alti¬ 
tude  where  A  -  arctangent  (AV/AZ)-  is  a  maximum*  as  shown  on  Fig.  1. 


Figure  1.  Sound  speed  profiles  defining  terms  for  the  Beta  parameter. 

(Left  figure  is  an  example  of  a  negative  Beta  parameter. 
Right  figure  depicts  a  positive  Bets  parameter.  Solid  lines 
are  the  sound  speed  vertical  profiles.  Dotted  line  is  the 
surface  sound  speed  value  for  comparison  with  other  levels. 
Small  dashed  lines  are  magnitudes  of  AV  and  AZ.  The  angles 
designated  A  are  the  maximum  values  of  arctangent  AV/AZ.) 

BIX 
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The  magnitude  of  the  Bata  parameter  along  a  particular  direction  la  computed 
from  the  maximum  valua  of  A  given  by  the  aquation: 


B  «  arctangent  13  *  (AV/AZ)  *  (R/C)]  (2) 

where 

B  "  Bata  paramatar  for  salactad  direction  (dag) 

R  ■  Dlatance  to  location  of  intaraat  (km) 

C  ■  3ound  apaad  at  the  aurfaca  (m/a) 

AV  *  Sound  apaad  diffaranca  related  to  A 

In  Fig.  1  (m/s) 

AZ  ■  Height  of  A V  above  ground  level  (km) 


Lorens'  alrblast  measurements  represent  the  positive  overpressure 
amplitude  or  the  alrblast  on  the  assumption  that  the  measurements  are  suffi¬ 
ciently  far-field  that  the  peak  positive  and  negative  overpressures  are 
approximately  equal.  The  weather  conditions  affecting  alrblast  propagation 
showed  a  significant  variation  In  the  vertical  temperature  gradient  and  wind 
velocity.  Resultant  Beta  parameters  ranged  from  a  +81  (strong  enhancement) 
to  a  -26  (moderate  attenuation).  A  least  squares  fit  to  the  data  has  the 
following  equation: 

L  *  103.1  +  B/5.3  (3) 

where 

L  ■  Instantaneous  peak  overpressure  (dB) 

B  -  Beta  parameter  as  previously  defined  (deg) 

Equation  3  has  a  standard  deviation  of  7.6.  Lorens  applied  the  results  from 
previous  detonations  with  yields  ranging  from  AS  to  540,000  kg  at  distances 
from  5  to  50  km  to  develop  the  yield  and  range  scaled  predictive  equation: 

L  -  103.1  +  B/5.3  +  20*L0G  [(S/i013)O,556*(W/ii0)0,4A4*(25/R)1,333]  (4) 
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where : 


L  ■  Iiaximum  peak  overpressure  (dB) 

B  »  Weather  parameter,  as  previously  defined  (deg) 

S  •  Surface  atmospheric  pressure  (rabar) 

W  *  TNT  equivalent  explosive  weight  (kg) 

R  *  Distancr  from  explosion  (km) 


Equation  4  represents  the  mean  expected  overpressure.  The  one  standard 
deviation  above  the  mean  is  obtained  by  adding  8  dB;  one  standard  deviation 
below  the  mean  is  obtained  by  subtracting  10  dB.  Lorenz'  model  was  developed 
with  measurements  in  units  of  decibels.  Decibels  are  an  arbitrary  scale  measure 
used  in  the  mechanical  measurement  of  sound.  This  arbitrary  scale  is  a  function 
of  the  ratio  of  the  instantaneous  overpressure  to  20  micropascals  which  is 
considered  the  threshold  of  human  hearing.  This  relationship  between  decibels 
and  pascals  can  be  seen  in  the  equation: 

L  -  20  log1Q  (p/Po) 

where : 

L  ■  Maximum  peak  overpressure  (dB) 

p  ■  instantaneous  overpressure  (Pa  x  10^) 

Pc  *  20  micropascals 

So  an  instantaneous  overpressure  reading  of  20  micropascals  is  equal  to  a 
maximum  peak  overpressure  of  zero  when  expressed  in  decibels.  The  current 
standard  of  airblast  intensity  is  units  of  pascals,  where  the  conversion  factor 
to  pascals  is  given  by: 

PK  -  0.00002  *  10  <L/20J  (Pa)  (5) 

Lorenz  stated  that  his  model  appeared  to  contain  the  core  of  a  fairly 
general  prediction  method.  However,  he  added  that  further  study  would  be  needed 
to  determine  the  applicability  >f  the  model  outside  the  range  of  data  from  which 
it  was  derived.  The  next  two  sections  demonstrate  how  well  the  model  performs, 
with  slight  modifications,  for  the  varied  AFWL  high- explosive  test  programs. 
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CONVENTIONAL  HIGH  EXPLOSIVES  BLAST  AND  SHOCK  (CHEBS)  SERIES 

This  section  describes  the  application  of  Lorens'  model. for  the  AFWL 
Conventional  High  Explosives  Blast  and  Shock  (CHEBS)  test  series.  The 
purpose  of  CHEBS,  conducted  near  the  southern  perimeter  of  Kirtland  APB,  was 
to  test  the  effects  if  conventional  Mark-83  1000-lb  bombs  on  protective  shel¬ 
ters  and  to  characterize  close-in  overpressures.  AFWL/WE  monitored  these 
tests  to  insure  that  unacceptable  airblast  would  not  extend  to  the  southern 
suburbs  of  Albuquerque,  approximately  15  km  from  the  test  site. 

For  all  CHEBS  shots,  the  airblast  Intensity  was  measured  at  a  location 
5  km  east-northeast  of  the  test  site.  Additional  measurements  at  7  and  10  km 
were  added  as  more  airblast  measuring  equipment  became  available.  Overpressures 
were  digitally  recorded  on  Valldyne  Differential  Pressure  Transducers,  Model 
P305D.  The  transducer  is  of  a  diaphragm  type  and  capable  of  measuring  in  a 
range  from  2  to  862  Pa,  with  the  manufacturer's  specifications  Indicating  an 
accuracy  of  ±  0.5  percent  at  full  scale;  i.e.,  ±  4.3  Pa.  Reference  8  gives 
a  complete  description  of  the  Valldyne  specifications  and  Fig.  2  shows  an 
example  of  the  recorded  airblast  waveform  from  the  Valldyne. 


Figure  2.  Example  of  Valldyne  plot  of  recorded  airblast  waveform. 
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Atmospheric  conditions  were  determined  using  the  regularly  scheduled 
ravinsonde  observations  from  the  Albuquerque  Airport,  approximately  10  km 
west  of  the  test  cite.  Test- time  vertical  temperature  profiles  were  estimated 
from  the  1200  Greenwich  .lean  Time  (GMT)  rawlnsonde  data,  adjusting  for  the 
observed  surface  temperature  at  shot-time,  normal  diurnal  variations,  and  any 
significant  advective  changes.  In  addition,  the  winds  aloft  were  measured  at 
shot  time  by  visually  tracking  a  10-g  pilot  balloon  (plbal)  with  a  theodolite. 
The  pibal  winds  could  not  be  determined  above  5000  ft  due  to  clouds  or  strong 
winds  limiting  the  elevation  angle  toward  the  Manzano  Mountains  to  the  east. 

For  winds  above  the  limit  of  plbal  data9  measurements  of  the  latest  available 
Albuquerque  rawlnsonde  winds  were  used  with  the  necessary  advective  adjustments. 

The  Mark-83  bombs,  which  were  placed  either  nose  into  the  ground  or  on 
the  side,  effectively  acted  as  a  point  surface  detonation  for  far-field  air- 
blast  intensity  purposes.  For  a  surface  burst,  the  blast  propagates  hemis- 
pherically,  rather  than  spherically  as  in  the  case  of  an  elevated  burst. 

Assuming  the  ground  surface  as  a  perfect  reflector,  the  resulting  overpressures 
are  approximately  equivalent  to  those  from  double  the  explosive  yield.  Thus, 
a  charge  weight  of  550  kg  TNT  equivalent  ( 2*275  kg)  was  used  as  the  explosive 
source  strength  for  the  Mark-83  bombs. 

The  results  of  the  CHEBS  series  are  listed  in  Table  2.  Weather  conditions, 
listed  in  Appendix  A,  resulted  in  the  Beta  parameter  ranging  from  -4  to  +11. 

All  overpressure  measurements  were  well  within  the  upper  bound  of  Lorenz* 
prediction  equation.  However,  the  mean  prediction  tended  to  be  low,  espe¬ 
cially  when  stronger  winds  aloft  resulted  in  a  positive  Beta  parameter. 
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TABLE  2.  CHEBS  RESULTS :  MEAN  AND  MAXIMUM  PREDICTIONS  USING  EQUATION  4 


Measurement 

Beta 

Beta 

Measured 

Mean 

Maximum 

Location 

Parameter 

Height 

Overpressure 

Prediction 

Prediction 

(Azlmnth/Range) 

(deg) 

(ft) 

(Pa) 

(Pa) 

(Pa) 

(deg/ka) 


7 

070/5 

+9 

2000 

81 

57 

145 

8 

070/5 

+11 

5000 

99 

60 

151 

9 

c/0/5 

-4 

1000 

44 

43 

109 

10 

070/5 

-2 

9000 

44 

45 

114 

12 

070/5 

-3 

7000 

45 

44 

112 

13 

070/5 

0 

1000 

71 

47 

119 

13 

070/10.4 

0 

1000 

30 

18 

46 

15 

070/5 

+3 

5000 

85  ‘ 

50 

132_ 

15 

070/10.4 

+8 

5000 

20 

22 

55 

15 

055/7.2 

+4 

5000 

30 

37 

82 

Lorenz'  mean  prediction  equation  was  adjusted  to  fit  the  data  and  given  by 
the  equation: 

L  -  106  +  B/5.3  +  20*LOG[(P/1013)°*556*(W/110)0*444*(25/R)1,333]  (6) 

where 


L  «  Peak  overpressure  (dB) 

B  *  Weather  parameter  (deg) 

P  ■  Ambient  pressure  (mbar) 

W  ■  TNT  equivalent  explosive  weight  (kg) 

R  >  Range  of  interest  (kn) 


Equation  5  is  used  to  convert  the  overpressure  to  Pascals.  The  maximum  over¬ 
pressure,  consistent  with  Lorens'  upper  bound*  is  obtained  by  adding  5  dB  to 
the  mean  predictive  value  of  Equation  6. 
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The  predictions  using  the  Equation  6  instead  of  Equation  4  are  listed 
in  Table  3.  Figure  3  depicts  the  linear  fit  of  the  predicted  overpressures 
to  the  observed  values  along  with  95  percent  confidence  intervals  (Ref.  9) 
for  a  future  single  response.  The  least  squares  fit  to  the  data  has  the 
equation:  Observed  »  -15.3012  -4-  (1.2208  *  predicted)  with  a  correlation 
coefficient  of  0.8801.  The  magnitude  of  the  confidence  intervals  is  wide  but 
reflects  the  small  sarple  size  (10  measurements)  rather  than  the  goodness  of 
fit.  Since  weather  conditions  were  partially  estimated  from  rawinsonde  data 
nearly  6  h  old,  the  predictions  using  Equation  6  in  the  CHEBS  series  were 
quite  reasonable.  The  next  section  shows  how  the  BOOM  prediction  technique 
performs  much  better  when  current  weather  data  are  available. 


TABLE  3.  CEEBS  RESULTS:  MEAN  AND  MAXIMUM  PREDICTIONS  USING  EQUATION  6  - 


Event  Measurement 

Beta 

Beta 

Measured 

Mean 

Maximum 

Location 

Parameter 

Height 

Overpressure 

Prediction 

Prediction 

( Az imuth / Range ) 

(deg) 

(ft) 

(Pa) 

(Pa) 

(Pa) 

(deg /km) 


7 

070/5 

+9 

2000 

81 

79 

145 

8 

070/5 

+11 

5000 

99 

82 

151 

9 

070/5 

-4 

1000 

44 

59 

109 

10 

070/5 

-2 

9000 

44 

61 

114 

12 

070/5 

-3 

7000 

45 

63 

112 

13 

070/5 

0 

1000 

71 

64 

119 

13 

070/10.4 

0 

1000 

30 

24 

46 

15 

070/5 

+3 

5000 

85 

70 

132 

15 

070/10.4 

+8 

5000 

20 

29 

55 

15 

070/7.2 

+4 

5000 

30 

44 

82 

11 


Figure  3.  Correlation  of  predicted  versus  measured  CHEBS  overpressures. • (The 
solid  line  is  the  linear  least-squares  fit:  Measured  *  -15.3012 
+  [1.2208  *  Predicted]  with  a  correlation  coefficient  of  0.8801. 
Dotted  lines  are  95  percent  confidence  intervals  for  a  future 
response. ) 
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I  CBM  SILO  SUPERHARE  MING  TECHNOLOGY  (ISST)  RESULTS 

The  purpose  of  the  Intercontinental  Beilis tic  Missile  Silo  Superhardening 
Technology  (ISST)  program,  conducted  at  e  renote  site  45  kn  southeast  of  Yuma, 
Arisons ,  vss  to  test  the  effects  of  simulated  nuclear  detonations  on  scale- 
eised  nlssile  silos.  To  attain  the  overpressures  necessary  to  sinulate  a 
nuclear  burst ,  the  explosives  (Irani te  60,  TNT  equivalent  0.90)  vere  laid  in 
a  mostly  circular  form  with  an  overburden  of  desert  alluvial  soil  placed  upon 
the  explosives.  The  configuration  was  designed  to  confine  most  of  the  explo¬ 
sives  effect  to  the  test-bed  area.  However,  for  airblast  prediction  purposes, 
the  design  led  to  the  problem  of  determining  the  portion  of  energy  emitted  to 
the  atmosphere.  This  section  describes  the  method  of  predicting  the  airblast 
attenuation  caused  by  an  overburden  and  presents  the  results  of  the  BOOM  model 
applied  to  airblast  predictions  for  the  ISST  series. 

Weathur  conditions  (wind  end  temperature  aloft)  were  gathered  from  ravin- 
sonde  data  measured  at  the  Yuma  Proving  Ground  (YPG) ,  40  km  northwest  of  the 
ISST  site.  Although  some  uncertainty  was  introduced  into  the  weather  data 
due  to  the  distance  from  the  test  site,  the  YPG  rewineonde  was  the  closest 
available  data.  The  data  were  always  recorded  within  30  min  of  test  time; 
therefore,  uncertainty  due  to  distance  was  minimised  by  the  timeliness  of 
the  Information.  This  problem  was  eliminated  for  the  Large  Size  #1  (LSI) 
test  (so  far,  the  largest  in  the  series)  when  the  YPG  weather  teem  recorded 
the  ravins onde  data  at  the  ISST  cite.  For  the  larger  ISST  tests,  weather 
people  from  the  Yuma  Marine  Corps  Air  Station  took  plbal  wind  observations 
at  the  site.  Since  temperature  is  e  much  more  conservative  meteorological 
quantity  than  wind  velocity,  the  combination  of  the  YPG  upper  air  temperature 
date  and  the  winds  aloft  observed  at  the  site  was  the  best  compromise  to  ade¬ 
quately  describe  the  state  of  the  atmosphere  for  airblast  propagation. 


Airblast  wu  again  recorded  on  the  AFWL  Validynes  for 
For  the  larger  tests,  additional  measurements  were  made  by 


the  ISST  events. 
Sandia  National 
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Laboratory  (SNL)  people,  who  aupportad  thaaa  taata  with  airblast  predictions 
and  measurements  (Refs.  10.  11,  and  12).  Table  A  lists  the  test  information, 
points  of  measurement ,  weather  parameters,  ana  the  measured  overpressures. 


TABLE  A.  ISST  TEST  RESULTS 


Example 


Location  Beta 

(Azimuth /Range)  Parameter 

_ _ <d«R) 


Measured 
Overpressure 
_ &*2 _ 


1 

355/17.1 

2 

355/3.6 

3 

360/13.8 

A 

3A1/3.8 

5 

359/19.3 

6 

340/14.4 

7 

355/16 

8 

355/3.3 

9 

337/11.3 

10 

340/9.6 

11 

337/5.3 

12 

337/15.6 

13 

330/9.7 

1A 

3A0/3.5 

15 

330/16.2 

16 

3A0/2.5 

17 

330/11.3 

18 

3A0/1A.A 

19 

300/2. 1 

20 

360/14.2 

21 

336/9.6 

22 

360/18.6 

23 

330/11.3 

2A 

336/9.7 

25 

3A0/9.6 

26 

355/ 1A. 5 

27 

330/9.7 

28 

3A0/3.5 

29 

3A0/1A.A 

30 

3A0/14.A 

31 

355/13.9 

32 

338/5.9 

33 

360/13.8 

-15 

2A 

-7 

133* 

-26 

11* 

-14 

345* 

-53 

28* 

-21 

1A 

+4 

40 

+1 

414* 

-10 

21 

-16 

20 

-5 

96 

-A6 

32 

-1A 

37* 

-7 

80 

+10 

58 

-1 

196 

+7 

38 

-36 

* 

-2 

373 

+3 

28* 

-33 

82 

-3A 

7* 

-16 

36 

-33 

75* 

-1A 

28 

-28 

14 

+6 

62* 

-5 

114 

-2A 

8 

-1 

6 

-AA 

39* 

+4 

190* 

-44 

49 

*  SNL  measurement. 
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To  determine  the  airblasc  attenuation  effect  of  the  ISST  overburden 
configuration ,  the  BOOM  predictive  equation  (Equation  6)  wee  eolved  for  the 
amount  of  exploaive  weight  required  to  produce  the  obeerved  overpressure  for 
the  range  of  interest  and  the  weather  conditions.  In  effect ,  thia  procedure 
eliminates  the  overburden  effect  by  treating  the  explosion  aa  a  point  source 
surface  detonation.  Therefore,  the  percentage  ratio  of  the  "apparent  weight" 
to  the  total  explosive  weight  describes  the  portion  of  energy  remaining  as 
airblast  after  the  explosive  energy  has  been  expended  to  remove  the  overburden. 
The  last  column  of  Table  5  shows  the  percentage  of  the  total  explosive  weight 
that  resulted  in  the  measured  overpressure  at  each  location.  For  the  series, 
values  ranged  from  less  than  1  percent  to  nearly  11  percent. 

J.W.  Reed  (Ref.  13)  reasoned  that  die  airblast  emanating  from  beneath  an 
overburden  should  depend  on  the  amount  of  explosives,  the  area  over  which  it 
is  laid,  the  overburden  mass,  and  the  time  to  venting.  However,  he  added  that 
the  time  and  character  of  the  exploaive  venting,  i.e.,  the  failure  of  he 
overburden  containment,  should  depend  on  the  texture  of  the  overburden  mate¬ 
rial.  Any  physical  assessment  of  that  characteristic  would  be  difficult. 

On  the  other  hand,  Reed  observed  that  tha  fractional  amount  of  energy  expended 
In  lifting  the  overburden  to  its  venting  point  depends  on  the  overburden  mass 
per  unit  explosive  area.  Thus,  the  vented  airblast  should  have  some  relation¬ 
ship  to  the  overburden  mass  per  unit  explosive  mass. 

As  prsviously  discussed,  the  purpose  of  the  ISST  overburden  was  to  induce 
large  overpressures  by  confining  most  of  the  explosive  energy  to  the  test-bed 
area.  Aa  depicted  in  Fig.  4,  the  explosives  we;:^  laid  in  a  nearly  circular 
area  with  the  overburden  built  to  a  prescribed  depth  over  the  explosives.  The 
depth  of  the  explosives  was  negligible  compared  to  the  overburden  height. 

The  edge  of  the  overburden  extended  beyond  the  explosives  extremity  to  a  dis¬ 
tance  equal  to  the  overburden  height  and  then  extended  to  the  ground  on  a 
2:1  grade.  This  configuration  was  designed  to  attain  an  equal  amount  of  over¬ 
burden  in  all  directions  from  the  explosives. 


15 


AFWL-TR-3.V-150 


TABLE  5.  ISST  APPARENT  WEIGHT  (X  OF  EXPLOSIVE  WEIGHT) 


Example 


Measurement 
Location  (lac) 


Apparent  Weight 
(X  of  Explosive  WelchtJ 


1 

4. 

17.1 

10.915 

2 

3.6 

2.707 

3 

13.8 

1.849 

4 

3.8 

2.220 

5 

19.3 

6.497 

6 

14.4 

4.999 

7 

16.0 

6.438 

8 

3.3 

6.347 

9 

11.3 

1.824 

10 

9.6 

4.477 

11 

5.3 

4.511 

12 

15.6 

3.627 

13 

9.7 

4.167 

14 

3.3 

3.163 

15 

16.2 

10.955 

16 

2.5 

4.734 

17 

11.3 

1.299 

18 

14.4 

3.012 

19 

2.1 

5.140 

20 

14.2 

4.707 

21 

9.6 

3.544 

22 

18.6 

2.633 

23 

11.3 

6.866 

24 

9.7 

3.140 

25 

9.6 

5.006 

26 

14.5 

3.112 

27 

9.7 

2.357 

28 

3.5 

3.681 

29 

14.4 

2. -841 

30 

14.4 

0.583 

31 

13.9 

3.458 

32 

5.9 

*>.285 

33 

13  8 

5.388 

>\ 


?  >  ■ 
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Figure  4.  Test-bed  configuration  for  the  ISST  series.  (Explosives  and 

overburden  ere  laid  In  a  nearly  circular  form.  The  dimensions  r, 
r  +  h,  and  r  +  3h  are  radii.  The  depth  of  explosives  Is  negli¬ 
gible  compared  to  the  overburden  height.) 


Many  approaches  were  tried  to  empirically  describe  the  airblast  attenua¬ 
tion  as  a  result  of  the  overburden.  The  best  relationship  resulted  when  the 
ratio  of  the  explosive  weight  and  the  total  overburden  mass  was  correlated 
with  the  mear  apparent  explosive  weight  for  each  individual  test.  The  entire 
overburden  mass  was  used  to  calculate  the  overburden  factor  with  the  assump¬ 
tion  of  a  haaispherically  radiating  blast  wave.  This  accounts  for  the  portion 
of  blast  energy  needed  to  move  the  intervening  overburden  in  all  directions 
before  venting  to  tbe  atmosphere.  The  overburden  factor  calculated  for  the 
ISST  series  is: 


where 


Z  -  W/(0*D) 


(7) 


Z  ■  Overburden  factor 
U  -  TNT  equivalent  explosive  weight 
0  "  Overburden  volume 
0  ■  Density  of  overburden  soil 


(kg  explosive/kg  overburden) 
(kg,) 

(n3) 

(kg/m3) 
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Tha  density  of  tha  Alluvial  soil  vaa  measured  aftar  aach  taat  and  found  to  be 
1845  hg/*"*,  i  5  parcant.  (Rote:  Equation  7  has  boon  taatad  for  only  thia 
density.) 

Tha  aquation  darivad  to  account  for  tha  ovarburdan  affacta  haa  tha 

fora: 


ZZ  -  -0.0146  ♦  (11.5157*Z)  (8) 

where 

ZZ  -  Effactlva  axploslva  weight  factor 
Z  -  Ovarburdan  factor  froa  Equation  7 

Tha  corralation  coafficiant  of  tha  linaar  fit  la  0.7660.  Figura  5  graphically 
shows  tha  laaat  aquaras  fit  of  tha  data.  Seattar  la  avid ant  in  tha  data; 
however,  tha  ovarburdan  factor  Z  doas  hava  a  poaitlva  linaar  corralation  with 
tha  apparant  explosive  vaight.  Aa  an  example,  if  an  aqual  ovarburdan  configu¬ 
ration  vaa  us  ad  vith  an  incraasad  axploalva  vaight*  Equation  8  corraetly  pre¬ 
dicts  an  incraasad  apparant  vaight  resulting  in  incraasad  airblast  affect. 
Conversely*  an  incraasad  ovarburdan  volume  vith  tha  case  axploslva  vaight 
vould  lead  to  a  decreased  airblast. 

The  ovarburdan  attenuation  (Equation  8)  vas  used  in  tha  BOOM  prediction 
aquation  (Equation  6)  to  calculate  tha  axploslva  vaight  for  tha  ISST  events. 
Tha  prediction  results  era  listed  in  Table  6.  The  95  parcant  confidence 
Halts  (Ref.  9)  vara  coaputad  for  a  single  future  response.  Based  on  tha 
results  of  33  overpressure  measurements ,  tha  upper  bound  for  predictions  in 
tha  ISST  series  is: 

FM  -  FK  4  43. 2638* [1.0303  ♦  ((FK-88.1212)2/389009. 5152] 1/2  (9) 

vhera 

FM  ■  Upper  bound  on  tha  BOOM  prediction  (Fa) 

PK  •  BOOM  naan  prediction  (Fa) 
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TABLE  6.  ISST  BOOM  OVERPRESSURE  PREDICTIONS 


Measured  Predicted 


Rxanple 

Maasur«ant 
Distance  (to) 

Bata 

(dea) 

Overpressure 

(Pe) 

Overpressure 

(P*> 

1 

17, 1 

-15 

24 

17 

2 

3.6 

-7 

133* 

152 

3 

13.8 

-26 

11* 

17 

4 

3.8 

-14 

345* 

444 

5 

19.3 

-53 

28* 

22 

6 

14.4 

-21 

14 

15 

7 

16.0 

+4 

40 

48 

8 

3.3 

+2 

414* 

372 

9 

11.3 

-10 

21 

33 

10 

9.6 

-16 

20 

17‘ 

11 

5.3 

-5 

96 

99 

12 

15.6 

-46 

32 

33 

-J 

9.7 

-14 

37* 

35 

14 

3.5 

-7 

80 

79 

15 

16.2 

+10 

58 

55 

16 

2.5 

-1 

196 

190 

17 

11.3 

+7 

38 

82 

18 

14.4 

-36 

7 

8 

19 

2.1 

-2 

373 

363 

20 

14.2 

+3  _ 

28a 

56 

21 

9.6 

-33 

82 

86 

22 

18.6 

-34 

7a 

9 

23 

11.3 

-16 

36 

27 

24 

9.7 

-35 

75* 

78 

25 

9.6 

-14 

28 

25 

26 

14,5 

-28 

14 

15 

27 

9.7 

+6 

62* 

99 

28 

3.5 

-5 

114 

138 

29 

14.4 

-24 

8 

8 

30 

14.4 

-1 

6 

11 

31 

13.9 

-43 

39a 

42 

32 

5.9 

+4 

190* 

183 

33 

13.8 

-43 

49 

42 

*  SNL  BUkuraatnta 
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As  shown  on  Fig.  6,  the  BOOM  prediction  technique  produced  excellent 
results,  i.e.,  a  0.9471  slope  of  the  linear  least-squares  fit  and  a  correla¬ 
tion  coefficient  of  0.9788.  Accurate  predictions  were  made  over  a  wide  range 
of  both  test  parameters  and  atmospheric  conditions.  Overpressure  predictions 
were  accurate  at  distances  from  2.1  to  18.6  km  at  overpressure  levels  ranging 
from  6  to  414  Pa.  Additionally,  weather  conditions  were  quite  variable. 
Vertical  gradients  of  both  temperature  and  wind  velocity  resulted  in  the 
Beta  parameter  varying  from  +10  to  -53.  However,  the  1SST  overburden  factor 
[explosive  weight/ (overburden  volume  *  overburden  density)]  was  only  between 
0.003  and  0.007.  Further  verification  is  needed  for  airblast  predictions 
using  overburden  factors  beyond  the  ISST  limits. 

A  computer  program  to  incorporate  the  BOOM  technique  was  written  in 
BASIC  programming  language  for  use  on  a  Radio  Shack  PC-2  portable  micro-  - 
computer.  See  Appendix  B  for  details  of  the  program. 
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2SST  PREDICTED  VS.  MEASURED  OVERP 


Correlation  of  predicted  and  measured  ISST  overpressures.  (The  solid 
line  is  the  linear  least-squares  fit:  Measured  »  -1.8014  +  (0.9471  * 
Predicted)  with  a  correlation  coefficient  of  0.9788.  Dotted  lines  art 
95  percent  confidence  intervals  for  a  future  response.) 
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CONCLUSION 


An  airblast  prediction  technique*  the  BOOM,  has  been  presented.  The 
BOOM  was  implemented  on  a  Radio  Shack  PC-2  portable  microcomputer  and  is 
particularly  applicable  for  airblast  predictions  at  remote  locations  where 
access  to  a  mainframe  computer  is  not  available.  The  BOOM  is  a  modification 
of  a  technique  developed  by  the  Naval  Surface  Weapons  Center  for  predicting 
airblast  intensities  from  500- lb  bomb  explosions.  The  refractive  effect  of 
the  atmosphere  is  incorporated  into  a  single  function,  as  opposed  to  ray¬ 
tracing  techniques,  to  determine  whether  airblast  intensity  will -be  amplified 
or  diminished  by  the  atmosphere.  Results  of  airblast  intensity  measurements 
from  1000-lb  bomb  explosions  demonstrate  the  validity  of  the  original  tech¬ 
nique.  An  overburden  factor,  based  upon  the  volume  of  soil  placed  upon  the 
explosives,  has  been  developed  to  predict  the  airblast  emanating  from  beneath 
an  alluvial  soil  overburden.  Airblast  predictions  using  the  BOOM  were  vali¬ 
dated  with  measurements  from  the  ISST  program.  The  overburden  airblast 
attenuation  factor  may  not  be  valid  for  ratios  of  explosive  weight  to  over¬ 
burden  volume  exceeding  the  range  of  the  ISST  series.  However,  the  general 
technique  should  be  useful  in  developing  attenuation  factors  for  explosions 
with  overburdens. 
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APPENDIX  A 

WEATHER  DATA  AND  SOUND  SPEED  PROFILES-CHEBS  SERIES 

This  appendix  lists  ths  atmospheric  conditions  and  tha  resultant  tabular 
and  graphical  vertical  sound  speed  profiles  for  the  CHEBS  high-explosive 
test  series.  The  listings  are  copies  of  print-outs  from  the  microcomputer 
-  used  to  run  the  program.  For  ease  of  readability  of  the  sound  speed  plots, 
an  enlarged  copy  of  the  first  plot  is  presented.  See  Section  4  for  the 
methods  of  obtaining  the  meteorological  data  and  the  assumptions  considered. 
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CHEBS7  11APR84 
*UEATHER  0A7A* 
HGT  TEnP  WIND 
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-6. 
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DIRECTION  70 

tSOUNO  SPEED  DATA* 


HGT  SOUND  DELTA 
SPEED  SPEED 

<kft)  «vs>  <rvs> 


0.0 

351.5 

1.0 

352.5 

1.0 

2.0 

353.7 

2.2 

3.0 

350.4 

-1.0 

4.0 

348.6 

-2.8 

5.0 

346.0 

-5.4 

6.0 

346.6 

-4.8 

7.0 

348.3 

-3. 1 

8.0 

352.6 

1.1 

9.0 

352.2 

0.7 

10.0 

352.3 

0.8 

Figure  A-l.  CHEBS  7  weather  data  and  sound  speed  profiles  toward  070  deg. 
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Figure  A-2.  Enlarged  example  of 
a  sound  speed  plot. 
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CMEBS8  3MAY84 
WEATHER  DATA* 


HGT 

TEnp 

UINO 

<KFT) 

<c> 

<KTS) 

0 

15 

270  10 

1 

12.6 

270  10 

2 

10.2 

275  16 

3 

10 

281  13 

4 

10 

280  25 

5 

3.3 

280  35 

6 

5.6 

230  35 

7 

1.3 

235  33 

8 

-1.6 

235  40 

3 

-5. 1 

300  41 

10 

-6.4 

235  43 

DIRECTION  70 


♦SOUND 

i  SPEED 

DATA* 

HGT 

SOUND 

DELTA 

SPEED 

SPEED 

<KFT) 

CfVS) 

cn/S) 

0.0 

345. 1 

1.0 

343.7 

-1.4 

2.0 

344.8 

-0.2 

3.0 

343.0 

-2.0 

4.0 

348.3 

3.1 

5.0 

352.2 

7.0 

6.0 

348.2 

3. 1 

7,0 

346.4 

1.7 

8.0 

344.6 

-0.5 

3.0 

341.5 

-3.6 

10.0 

342.7 

-2.3 

Figure  A-3.  CHEBS  8  weather  data  and  sound  speed  profile  toward 
070  deg. 
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CHE339  20AUG34 

♦LEATHER  DATA* 

HuT  TErP  LINO 
<KFT>  (C>  <KTS> 


0  22  219  4 

1  19.5  240  5 

2  12.1  240  5 

3  14.9  230  2 

4  12.6  2-10  6 

5  13.3  220  11 

6  8.2  230  11 

2  6.1  210  12 

3  4  213  14 

9  1.9  213  19 

10  -0. 1  210  21 


DIRECTION  20 


♦SOUND  SPEED  DATA* 

HST  SOUND  DELTA 
SPEED  SPEED 
(KFT >  <rvs>  <rvs> 


0.9  346.1 

1.3  345.5  -0.5 

2.0  344.1  -1.9 

3.0  343.2  -2.4 

4.0  342.0  -4.1 

5.0  342.4  -3.6 

6.0  341.6  -4.5 

2.0  339.2  -6.3 

8.0  339.2  -6.8 

9.0  339.9  -5.1 

10.0  339.4  -6.6 


Figure  A-4.  CHEBS  9  weather  data  and  sound  speed  profile  toward 
070  deg. 
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CHESS 10  27SEP84 

WEATHER  DATA* 

HGT  TEnP  UIND 
CKFT )  <C)  (KTS) 


0  12  330  5 

1  10.7  330  4 

2  9.4  340  4 

3  8.1  335  9 

4  7  325  11 

5  4.2  320  11 

6  1.7  310  12 

7  0  305  12 

8  -1  305  13 

9  -1.9  285  15 

10  -3.5  285  15 


DIRECTION  70 


*SOUNQ  SPEED  DATA* 

HGT  SOUND  DELTA 
SPEED  SPEED 
<KFT )  crvs>  <n/S> 


O. 0  339.1 

1.0  338.2  -0.8 

2.0  337.1  -1.9 

3.0  336.7  -2.3 

4.0  337.1  -1.9 

5.0  335.9  -3.2 

6.0  335.5  -3.6 

7.0  334.9  -4.1 

P. ff  334.6  -4.5 

9.0  336.4  -2.6 

10.0  335.5  -3.6 


Figure  A-5.  CHEBS  10  weaCher  data  and  sound  speed  profile  toward 
070  deg. 
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CFEBS12  14FEB85 
LEATHER  DATA* 
HGT  TEMP  UIND 


<KFT) 

<C> 

<KTS) 

0 

7 

360 

5 

1 

4.3 

360 

5 

2 

0.8 

340 

7 

3 

-1.5 

330 

3 

4 

-1.9 

315 

8 

5 

-3.  1 

315 

13 

6 

-4.4 

320 

18 

7 

-5.7 

320 

24 

8 

-7.4 

330 

28 

3 

-3.  1 

340 

32 

10 

-10.3 

i  340 

35 

DIRECTION  70 


tSOUNO 

I  SPEED 

DATA* 

HGT 

SOUND 

DELTA 

SPEED 

SPEED 

CKFT7 

<rvs) 

<f1/S> 

0.0 

334.8 

1.0 

333.2 

-1.6 

2.0 

331.3 

-2.8 

3.0 

331.3 

-3.5 

4.0 

332.0 

-2.8 

5.0 

332.3 

-2.5 

6.0 

331.8 

-2.3 

7.0 

332.1 

-2.7 

8.0 

323.4 

-5.4 

3.0 

325.3 

-8.3 

10.0 

324.8 

-10.0 

Figure  A-6 .  CHEBS  12  weather  data  and  sound  speed  profile  toward  070 
deg. 
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DIRECTION  70 


CHESS 13  28FESSS 


♦SOUND  SPEED  DATA* 


♦WEATHER  DATA* 


HGT 

<KFT 

TEMP 
>  <C) 

WIND 

CKTS) 

HGT 

(KFT) 

SOUND 

SPEED 

<0/S> 

DELTA 

SPEED 

(O/S) 

0 

11 

130 

5 

0.0 

339.3 

I 

8.7 

210 

6 

1.0 

339.0 

-0.3 

2 

5.3 

230 

6 

2.0 

337.8 

-1.4 

3 

3. 1 

280 

10 

3.0 

337.7 

-1.6 

4 

0.8 

300 

15 

4.0 

336.8 

-2.5 

5 

~1.3 

300 

13 

5.0 

334.8 

-4.4 

6 

-3.4 

280 

10 

6.0 

333.7 

-5.5 

7 

-5.5 

270 

10 

7.0 

332.8 

-6.5 

a 

-7.7 
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10 

8.0 

331.7 

-7.6 

9 

-9.9 
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12 

9.0 

331.4 

-7.9 

10 

-11.9 

240 

15 

10.0 

331.5 

-7.7 

Figure  A-7.  CHEBS  13  weather  data  and  sound  speed  profile  toward 
070  deg. 
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CHEBSI5  25APR85 
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DIRECTION  20 


tSOUNO  SPEED 

OAT  At 
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SOUND 
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SPEED 

SPEEO 

<KFT) 

CIVS) 

ervs) 

0.0 
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1.0 

345. 0 

-1.4 

2.0 

343.5 

-2.9 

3.0 

344.4 

-2.0 

4.0 

343.6 

-2.8 

5.0 

348.8 
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6.0 

348.0 
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347.7 

1.2 

8.0 

347.3 

0.8 

9.0 
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DIRECTION  55 


tSOUND  SPEED 

OATAt 

H6T 

SOUND 

DELTA 

SPEED 

SPEED 

CKFT) 

<M/S> 

<rvs> 

0.0 

347.8 

1.0 

346.5 

-1.2 

2.0 

345.3 

-2.4 

3.0 

347.3 

-0.4 

4.0 

346.7 

-1.0 

5.0 

349.7 

1.9 

6.0 

347.5 

-0.2 

7.0 

345.5 

-2.2 

8.0 

345.4 

-2.3 

9.0 

348.4 

-1.3 

10.0 

344.5 

-3.2 

Figure  A- 8.  CHESS  15  weather  data  and  sound  speed 
profiles  toward  070  deg  and  055  deg. 
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APPENDIX  B 

DESCRIPTION  OF  THE  BLAST  OPERATIONAL  OVERPRESSURE  MODEL  (BOOM) 

This  appendix  describes  the  AFWL/WB  BOOM  used  to  predict  far- fie Id 
airblest  intensity.  The  program  is  written  in  BASIC  programming  language 
for  use  on  a  Radio  Shack  TRS-80  PC-2  portable  microcomputer.  The  BOOM 
program  has  many  options.  Airblast  intensity  is  computwi  for  either  a 
standard  surface  detonation  or  a  buried  explosion.  The  model  uses  an 
attenuation  factor  to  predict  the  airblast  emitted  from  beneath  an  alluvial 
soil  overburden.  Weather  information,  to  compute  the  vertical  sound  speed 
profile,  may  be  input  in  the  following  trays:  (1)  rawinsonde  pressure  levels 
and  temperatures,  plus  winds  at  standard  1000  ft  intervals;  (2)  rawinsonde 
pressure  levels  and  temperature,  plus  pilot  balloon  azimuth  and  elevation- 
angles  to  compute  pibal  winds;  and  (3)  temperature  and  wind  at  specified 
heights.  The  output  statements  and  the  sound  speed  plot  subroutine  are 
written  for  the  PC-2  printer  having  only  an  18  character  field  width. 
Therefore,  users  of  the  program  will  probably  need  to  modify  these  statements 
for  use  on  other  machines.  On  the  following  pages  are  listings  of  the  pro¬ 
gram  variables  and  the  BOOM  program,  and  a  description  of  the  program  flow. 
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BOOM  VARIABLE  LISTING 

-  trigonometric  functions  to  compute  tbs  U  coupon* a c  of  ths  plbsl  wind 
AV.BV  -  trigonometric  functions  to  cosputs  ths  V  cosponont  of  ths  plbsl  wind 
AZ  -  azimuth  of  ths  pilot  bslloon  st  ssch  tissd  measurement  (dsg) 

B  "  maximum  Bsts  parameter  for  s  givsn  dir  set  ion  sod  rang*  (dsg) 

BH  -  altituds  of  maximum  Bsta  parameter  (left) 

CH  -  ov  cion  to  print  vortical  sound  apssd  plot 
00  -  difference  batwaan  dlrsction  of  interest  and  wind  direction 
at  each  height  (dsg) 

DI  ~  directions  for  alrblast  predictions  (dag) 

ON  -  number  of  directions  for  alrblast  predictions 
DV  -  diffarsnes  batwaan  surface  sound  spaed  and  sound  speed 
at  each  height  (s/s) 

DZ  -  geopotential  thickness  between  each  rawlnsonde  level, 
based  on  the  hydrostatic  ai  <  oxiaatloo  (ft) 

EL  -  elevation  of  tha  pilot  balloon  at  each  tided  observation  (deg) 

F  -  linear  interpolation  factor  to  compute  tesperature  at 
the  height  of  each  wind  data  height 
H  -  height  of  sound  speed  calculations,  based  upon  input  height 
of  wind  and  temperature  data  (kft) 

HD  -  height  difference  between  successive  rawlnsonde  pressure  levels  (ft) 
II  -  counter  for  pibal  wind  calculations 
Id  -  pibal  observation  Interval  (s) 

IT  ~  total  time  elapead  since  pilot  balloon  release  (s) 

LA  -  previous  pilot  balloon  azimuth  ;used  to  compute  pibal  winds  (deg) 

LE  -  previous  pilot  balloon  elevation  ;  used  to  compute  pibal  winds  (deg) 
LL  -  number  of  weather  data  levels 

M  -  Beta  parameter  for  each  direction,  range  and  height  (deg) 

MM  -  number  of  pilot  balloon  azimuth  and  elevation  observations 
9H  -  height  of  pilot  balloon  at  each  observation  (ft) 

MR  -  number  of  ranges  for  airblast  prediction  along  each  direction 
P  *  rawlnsonde  pressure  level  (mbar) 

?0  -  pibal  wind  direction  (deg) 

PH  ~  previous  height  of  pilot  balloon;  used  to  calculate  pibal  winds  (ft) 
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PK  -  Man  ovarprassura  froa  BOOM  (Pa) 

PM  -  maximum  ovarpraaaura  froa  BOOM  (Pa) 

PS  -  plbal  wind  apaad  (kt) 

R  -  rangaa  for  airblaat  pradictlon  along  aach  dlraction  (ka) 

RH  -  gaopocantlal  height  of  aach  rawlnaonda  praaaura  level  (ft) 

RL  -  nuabar  of  rawlnaonda  praaaura  lavala 

RR  -  rlaa  rata  of  pilot  balloon  dapandlng  on  tlM  alnea  ralaaaa  (ft/a) 
S  -  atmospheric  aurfaca  praaaura  (abar) 

SZ  -  pilot  balloon  alaa  (10  or  30  g) 

T  -  caaperatura  at  aach  rawlnaonda  praaaura  laval  (°C) 

TD  -  caaperatura  dlffaranca  between  successive  rawlnsonde 
praaaura  levels  (°C) 

TE  -  option,  surface  o;  burled  explosion 

TK  -  rawlnsonde  temperature  converted  to  Kelvin  (K) 

TN$  -  beading,  any  desired  entry 

TP  -  option  for  wind  input,:  rawlnaonda  or  plbal 

TT  -  caaperatura  at  the  level  of  wind  data;  either  actual  or  computed 
froa  rawlnsonde  data  (°C) 

TW  -  option  for  type  of  weather  data  input 

U  -  east-west  plbal  wind  component;  used  to  compute  plbal  wind 
V  -  north-south  plbal  vied  component;  used  to  compute  plbal  wind 
VS  -  sound  speed  along  the  selected  dlraction  (m/s) 

W  -  TNT  equivalent  apparent  explosive  weight  (kg) 

HD  -  wind  direction  at  aach  haight  (dag) 

VS  -  wind  speed  at  aach  haight  (kt) 

VT  -  TNT  equivalent  explosive  weight  (kg) 

XH  -  height  of  layer  averaged  plbal  wind  (ft) 

Z  -  overburden  factor  (kg  exploslves/kg  soil  overburden) 

ZZ  «  effective  explosive  weight  factor  (resulting  from  overburden) 
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BOOM  OOMPUTBK  PtOGIAM  LISTING 


l : "BOOM" 

10:OIR  P<10), T( 10 
),  TKU0),  RH(  10 
),  WD(21  >,  USC21 
) 

1 1 : DIR  TT<21),US< 
3,  21),  0U<3,21) 
,  01(3),  RC3,  8), 
PK<3,  8),  NR<8) 
15: DIR  B<3,  8),PR< 
3,  8),H(2i),BH( 
3,  8 ) :  UA I T  100 
30: INPUT  "TEST  LO 
CATION,  DATE/ T I 
RE” 5  TN* 

32: INPUT  "EXPLOSI 
UE  HEIGHT (KG) " 
JUT 

34: INPUT  "SFC  PRE 
SSCRBV'jS 
36: INPUT  "EXPLOSI 
ON?  SFC- 1, BUR I 
ED-2" J  TE  • 

38: IF  TE*1THEN 
LET  U=2*UT: 
30TC  49 

40: INPUT  "OUERBUR 
DEN  FACTOR  ";Z 
42 :  2Z=- .  0 1 46-*-  (11. 
5157*Z):U*WT*Z 
Z 

49: INPUT  "UANT  PL 
OT?  YES-1  NO-2 

M  •  Q|-| 

50: INPUT  "#  OF  DI 
R  OF  I NT"; ON 
60: FOR  1  =  1  TO  DN 
70: INPUT  "01 RECTI 
ON" *01 <I> 

71  .“INPUT  "#  OF  RA 
N3ES  ";NRU) 

72 s  for  J*iro  NRU 
) 

73 i INPUT  " RANGE CK 

R)  “;R<i,j) 

75* NEXT  JiNCXT  I 


90: PRINT  "TYPE  UE 
ATHER  DATA?" 

92:  INPUT  !,RAUIN-1 
,  HGT,  TERP,  UNO- 
2" J TU 

94:  IF  T 14*2 THEN 
GOTO  1170 
100: INPUT  "#  OF  RA 
WIN  LEUELS  "JR 
L 

120: FOR  1*1  TO  RL 
130: INPUT  "PRE3SUR 
ECRB)  ";P(I ) 
140:  INPUT  "TERPCO 

" ;  t  c  i  ) 

156: TK<  I  >-*TC  1 >+273 

.2 

160: NEXT  I 
220: PR I NT  "WIND  DA 
TA  NEEDED" 

230: INPUT  "W*D  DAT 
A?  RAUIN-1JPI9 
A:  -2";  TP 
240: IF  TP* 1  THEN 
GOTO  900 

250: INPUT  "BALLON 
SIZECGR)  "JSZ 
26e: INPUT  ' OBS  TIR 
£  INTERUAL  "JI 
S 

270: INPUT  ’#  PIBAu 
OCS  "JHR 

280: IT*0:PH*0:LA*0 
*'.£*. 01745: Cl* 
8.33:C2=7.«7:C 
3*7:RAr.*Jl/!B0: 
11*0 

290: IF  SZ* 10 THEN 
GOTO  310 

300 :  C  l  -*C  U*2 :  C2-C2* 
2t.C3«C3*2 
c 10: LPRINT  "  PIB 
AL  DATA" 

315:LF  JULPRINT  "H 
GT  <FT  >  UiNO< 
KT>" 

340; L*  2 
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353:  17*1  7*IS:RR*C2 
: 11*11+1 

363: IF  I T<=75TH£N 
LET  RR*C1 
370: IF  IT>240THEN 
LET  RR*C3 

383 : NH*PH*RR# I S :  XH 
s*<NH+PH>/2 
390: INPUT  "AZIMUTH 
ANGLE  JAZ 
400: INPUT  "ELEUATI 
ON  ANGLE  " » EL 
402 : AZ=AZ*RAO : EL  *E 
L*RAO: RADIAN 
410: AU*(NH#SIN  AZ> 
/TAN  SL'.8U*(PH 
*SIN  LAJ/TAN  L 
E:U*CAU--3U)/IS 
420: AU=(NH*COS  AZ> 
/TAN  EL:0U=<PH 
*CGS  LA)/ TAN  L 
E:U=CAU-8U)/IS 
430:P3=CSQR  CU^Z+U 
^2)>/3,2808 
440: IF  U<=0T^en 
GOTO  460 

450: PD=ATN  CLVUVR 
AD+180: GOTO  47 
0 

460:PO=ATN  CU/U>/R 
AD 

470: IF  PD<«0THEN 
LET  PO»PO+360 
480:  LA*AZ:LE*El.:  PH 
*NH 

490: LPRINT  USING  " 
######" > XHJ 
500: TAB  9: LPRINT 
USING  "####"} P 
D: : TAB  13: 
LPRINT  USING  " 
#*ir}PS*2 
510:  IF  IKKMTHEN 
GOTO  350 
520: LF  4 

300: INPUT  "*  OF  UN 
D  LEUELS  "ILL 


910: FOR  1*1  TO  L^ 

320: INPUT  "HE IGHTC 
KFT)  ";h<d 
330: INPUT  "UNO  DIR 
"  ;UD( I ) 

S40: INPUT  "UNO  SPD 
(KT)  "JUSCI) 

950: US (I >=US( I )t.5 
960: NEXT  I 
1000: RHC 1 >-0 
1010: FOR  I»2TO  RL 
1030: DZ=-48. 0373* 
ctkcd+tkci- 
1 ) )#LN  (PCI) 
/PC I-l)> 

1 040 : RH  C I ) SRH  C I - 1 
)+DZ 

1050:  NEXT  I 
1060:C*0:NN*RL-1 
1080: FOR  1*1  TO  NN 
1090:HD”RHCI+1)-R 
HC I >: TD=T  C I  + 
D-TCI) 

11 10: FOR  J*1T0  LL 
1120: IF  CHCJ > #100 
0XRHCI  >THEN 
GOTO  1161 
1130: IF  CHCJ>#100 
0) >RHC I+l ) 
THEN  GOTO  11 
62 

1140: F*(H( J )#1000 
-RHC  I  >  )/’HD:  C 
»C+l : TT  CC)=T 
D*F+TCI> 

1161: NEXT  J 
1162:  NEXT  I 
1165:  GOTO  1200 
1170: PRINT  "ENTER 
HT,  TEMP,  UIN 
D" 

1172*  INPUT  "HOU  PI 
ANY  LEUELS?” 

*  J  LL 

1174: FOR  I*1T0  LL 
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1178: INPUT  "ENTER 
HEIGHT  (KFT) 
** ;  h<  i  > 

1178S INPUT  "ENTER 
TEHP<C>  ";T 
TCI) 

1188: INPUT  "ENTER 
UNO  OIR  " JU 
0(1) 

1182: INPUT  "ENTER 
UNO  SPDCKT) 
••JUSCI) 

1 1 84 : US  C I ) =US Cl)* 
.5 

1186: NEXT  I 
1200: FOR  K*1  TO  DN 
1210: DEGREE 
1250: FOR  0*1  TO  LL 
1250: 00*01 (K)-UOC 
O) 

1270:  USCK,  0)*331. 
5*SQR  Cl+TTC 
n)/,273. 2)-US 
cmscos  (00) 

1280:  DUCK,  f1)*USCX 

,  id-usck,  l) 

1230: NEXT  n 
1300: NEXT  K 
1310: LPRINT  TNt 
1315: LF  1 
1316: LPRINT  "  *W 
EATHER  DATA* 

M 

1317: LF  1 
1320: LPRINT  "  HGT 
TETIP  UIN 

0" 

1322: LPRINT  " (KFT 
)  CC)  CKTS 
)" 

1325:  LPRINT  " - 


1330: LF  1 

1340: FOR  N*1 TO  LL 


1 360: LPRINT  USING 
"44#. 4" » HCN) 

1361: TAB  6: LPRINT 
USING  "444.4 
" ;  TTCN) ; 

1362: TAB  11: 

LPRINT  USING 
"4444" ; UOCN) 

1363: TAB  15: 

LPRINT  USING 
"444" ;US(N)¥ 
2 

1380: NEXT  N 

1391 SLF  2 

1400: FOR  K*1 TO  ON 

1420: LPRINT  "DIRE 
CTION  "} 

1425: LPRINT  USING 
"4444"; 01 (K) 

1430: LF  1 

1434: LPRINT  "*SOU 
ND  SPEED  DAT 
A*" 

1435: LF  1 

1440: LPRINT  "  HGT 
SOUND  DEL 
TA" 

1441: LPRINT  " 

SPEED  SPE 
ED" 

1442: LPRINT  " (KFT 
)  (0/S)  (0/ 

S)" 

1443: LPRINT  " - 


1500: FOR  1*1  TO  LL 
1510SLPRINT  USING 
"444. 4" | HC I ) 

1530: TAB  6: LPRINT 
USING  "4444. 
4"|USCK,  I); 
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1530: TAB  13: 

lprint  using 

"444.  4"  ;QU(K 
,  I) 

1560: N£XT  I 
1930: FOR  1*1  TO  NR 
<K> 

1938: BHCK,  I>*0:B' 

K,  I  >*-1000 
1939:B(K,  I >*-1000 
1950: FOR  J*2T0  LL 
1960: DEGREE 
1970:M*ATN  C'(3*R 
(K,  I>*OU(K,  J 
>>/OJS(K, 1>* 
H(J>*.3048>> 

> 

1979: IF  n>BCK,  I> 
THEN  LET  BH< 
K,  I  >*H( J> 
1980:  IF  n>B(K*  I) 
THEN  LET  B(K 
,  i>*n 

1990: NEXT  J 
2000: L 1*1 05*B(K» I 
>/5.3+20*LOG 
C<CS^1013>^. 
556>*<<U/110 
>/'.444)*(<25 
✓R(K,  I >  >^1 . 3 
33  >  > 

2005: PK(K,  I >*.000 
02*10^<L1'20 
> 

2008: IF  TE*2THEN 
GOTO  2030 
2010: L2*l 1 1+B(K»  I 
>/5.3+20*LOG 
<  (  <S/'1013>^. 
556>*CCU/110 
>^.444>»(<25 
✓R<K,  I  >>*1.3 
33)  > 

2015: PriCK,  I >*.000 
02*10*<L2/20 
> 

2020: GOTO  2040 


2030: PH (K»  I >*PK(X 
,  I >*43. 2838* 
SQR  (1.0303* 
(C(PK(X, I >-3 
8. 1212>*2>'3 
89009. 5152>> 
2040: NEXT  I 
2060: LF  3 
2062: LPRINT  "  *¥ 

PREDICTIONS# 
#":LF  1 

2070: LPRINT  "RANG 
E  OUER-PRE 
SS" 

2072: LPRINT  "(Kfl) 
(PA>  " 

2073: LPRINT  " - 


— ":LF  1 

2080: FOR  1*1  TO  NR 
(K> 

2090: TAB  l: LPRINT 
USING  "444.4 
" ;  R(K>  I  >  j 
2110: TAB  9ILPRINT 
USING  "44444 
4.4";PK(K,  I> 
2115: TAB  6:LPR1NT 
"OAX"; 

21 16: TAB  9:LPRINT 
USING  "44444 
4.4"jPn(K,  I> 
2120: LF  1 
2122: LPRINT  "MAX 
BETA  "JBCK, 
I> 

2123: LPRINT  "HGT ( 
KFT)  "JBHCK, 
I> 

2124: LF  2 
2125: NEXT  1 
2126: IF  CH*1THEN 
GOSUB  3000 
2122: TEXT  :LF  3 
2130: NEXT  K 
2140: END 
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3000: ren  plot  ro 

UTINE 
3012: LF  3 
3820JLPRIN7  "*SOU 
NO  SPEEO  PLO 
T*” 

3030: LPRINT  H  DIR 
ECTION”;OI (K 
) 

3043: TEXT  :LF  23 
3080: CSI2E  1: 

LPRINT  ”-24 
-18  -12 
6  0  6 
12" 

3030: LF  1 : CSI2E  1 
4010:LPRINT  "OELT 
ft  SPEED  (M/’S) 
sUITH  RESPEC 
T  TO  SFC": LF 
-4: GRftPH 

4040:  LINE  (0,  0)-( 
180,  0),  8: 
CSIZE  It 
LPRINT  ”  0” 

4050:  LINE  (0,80)- 
(180,  90),  8: 
LPRINT  ”  4” 

4060: LINE  (0,160) 
-(180,  160),  8 
: LPRINT  "  8 

M 

4070:  L I  Nf.  .0,240) 
-(180,  240),  8 

: LPRINT  "  12 
•• 

4080: LINE  (0,320) 
-(180,  320),  8 
: LPRINT  "  16 

M 

4090: LINE  (0,400) 
- ( 1 80,  400),  8 
ILPR!  -  “ 


4100:  LINE  (130,40 
0)-(180,  0),  3 
:LINE  (150,0 
)-( 130,  400), 

a 

4110: LINE  (120, 40 
0)-(120»  0),  3 
: LINE  00,0) 
-(90,  400),  8 
4120: LINE  (60,400 
)-(6Q,  0),  8: 
LINE  (30,0)- 
(30,  400),  8 
4130: LINE  (0,400) 
-(0,  0),  9 

4130: FOR  1*1  TO  LL 
:H(I)*H(1)*2 
0: DU(K, I )*DU 
(K,  I )*5: NEXT 
I 

4160: GL CURSOR  (12 
0,  0) :  SORGN  : 
FOR  J*2T0  LL 
4170: LINE  (DU(K, J 
-1 )»  H( J-l ) )- 
OU(K,  J),  H(J 
)),  0 

4 130:  NEXT  J  , 
4183: GLCURSOR  (90 
,  150): ROTATE 
3: LPRINT  "ftL 
TITUDE  (KFT) 

M 

4190: FOR  I*1T0  LL 
:H(I)*H(I)/2 
0: OU(K,  I)*OU 
(K, I)/3:NEXT 
I 

4195: TEXT  : CSIZE 
2:LF  20 
3000:  RETURN 
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